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Abstract 

The spin network simulator model represents a bridge between (generalised) 
circuit schemes for standard quantum computation and approaches based on no- 
tions from Topological Quantum Field Theories (TQFTs). The key tool is pro- 
vided by the fiber space structure underlying the model which exhibits combina- 
torial properties closely related to SU{2) state sum models, widely employed in 
discretizing TQFTs and quantum gravity in low spacetime dimensions. 
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We discuss here some aspects of a novel setting for quantum computa- 
tion pP [2] based on the (re) coupling theory of SU{2) angular momenta 
On the one hand, the "spin network simulator" can be thought of as a 
non-Boolean generalisation of the quantum circuit model, with unitary gates 
expressed in terms of: 

i) recoupling coefficients (3nj symbols) between inequivalent binary coupling 
schemes of A = (n + 1) 5'f/(2)-angular momenta (j-gates); 

ii) Wigner rotations in the eigenspace of the total angular momentum (M- 
gates). The picture does contain the Boolean the particular case when 
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all N angular momenta are spin i. The combinatorial structure of the model, 
on the other hand, closely resembles SU{2) "state sum models" employed in 
discretized approaches to 2 and 3-dimensional Topological Quantum Field 
Theories |5j (TQFT) and Quantum Gravity ^ . An explicit mapping relating 
the spin network to the topological approach to quantum computation based 
on modular functors [Z] can be established as illustrated in Ref. 2, Section 
6. 

The basic ingredients of the kinematics of the simulator are computa- 
tional Hilbert spaces (to be defined below) equipped with all unitary gates of 
types i) and ii) we may consistently apply to transform vectors belonging to 
such spaces. A specific computation will be implemented by picking up an 
"input" state to get an "output" generated by acting with a particular (finite) 
sequence of elementary gates. An elementary operation will be carried out 
in one unit of time since it is assumed that the intrinsic time variable of the 
simulator has a discrete character. 

The architecture of the spin network is modelled as an SU{2) fiber space 
structure over a discrete base space V 

{V, su{2y\ (1) 

which encodes all possible computational Hilbert spaces as well as gates for 
any fixed number iV = n + 1 of incoming angular momenta. 

The base space V = {v{b)} represents the vertex set of a regular, 
2(n — l)-valent graph 0„(V,i?) with cardinality \V\ = {2n — 1)!!. There 
exists a 1:1 correspondence 

Mb)} ^ {Him (2) 

between the vertices of (5n{V,E) and the computational Hilbert spaces of 
the simulator. For a given value of n, 7i;[(b) is the simultaneous eigenspace 
of the squares of 2?7, + 1 Hermitean, mutually commuting angular momentum 
operators, namely 

Jl; J 2, Js;---; Jn+1 = {Jj}; Jl-|-J2 + J3+--- + Jn+1 = 

Ki, K2, K3, . . . , K„_i = {K4 (3) 

and of the operator (the projection of the total angular momentum J along 
the quantization axis). The associated quantum numbers are ji,j2, ■ ■ ■ ,jn+i 
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; J; ki,k2,..., kn-i and M, where —J < M < in integer steps. If 
W'^ (g) ■ ■ ■ (g) Ti-'^+i denotes the factorized Hilbert space, namely the 

{n + l)-fold tensor product of the individual eigenspaces of the (Jj)^ 's, the 
operators K/j's represent intermediate angular momenta generated, through 
Clebsch-Gordan series, whenever a pair of Jj's are (binarily) coupled. As an 
example, by coupling sequentially the Jj's according to the scheme (■ ■ ■ ((Ji + 

J2) + J3) H \- Jn+i) = J - which generates (Ji + J2) = Ki, (Ki + J3) = 

K2, and so on - we should get a binary bracketing structure of the type 
(■ ■ ■ iiiW' ® W^)ki ®W'')k2® ■ ■ ■ ® ^^"+Ofcn-i)j, where we add an overall 
bracket labeled by the quantum number of the total angular momentum J. 
Note that, as far as jVs quantum numbers are involved, any value belonging 
to {0,1/2,1,3/2,...} is allowed, while the ranges of the fc/j's are suitably 
constrained by Clebsch-Gordan decompositions {e.g. if (Ji + J2) = Ki ^ 
lii — i2\ < ki < ji + ^2)- We denote a binary coupled basis of {n + 1) 
angular momenta in the JM-representation and the corresponding Hilbert 
space appearing in (j21) according to 

{ I [Jl, J2, J3, • • . , Jn+l]' ; klkl..., kl_^ ■ JM), -J<M< J} 

= niib) = span{ |b;JM)„}, (4) 

where the string inside [ji, 72 , is, • • • , jn+i]^ is not necessarily an ordered one, 
b indicates the current binary bracketing structure and the khS are uniquely 
associated with the chain of pairwise couplings selected by b. For a given 
value of J each 7i;[(b) has dimension (2J + 1) over C and thus there exists 
one isomorphism 

niib) c2^+^ (5) 

for each admissible binary coupling scheme b of (n + 1) incoming spins. The 
vector space C^'^"'"^ is interpreted as the typical fiber attached to each vertex 
v{b) eV oi the fiber space structure (H)) through the isomorphism (0). 

For what concerns unitary operations acting on the computational Hilbert 
spaces OD), we examine first j-gates associated with recoupling coefficients 
(3nj symbols) of SU{2) as anticipated in i). It can be shown jS] that any such 
coefficient can be splitted into "elementary" j-gates, the Racah transforms, 
possibly apart from phases and weight factors. A Racah transform applied 
to a basis vector is defined formally as 

n : \...{{ab),c)f...;JM) ^ \ . . .{a{bc),)f . . JM), (6) 
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where we are using here Latin letters a,b,c, . . . to denote both incoming (jj 's 
in the previous notation) and intermediate {kh 's) spin quantum numbers. 
The exphcit expression of (jHl) reads 

|(a(6c)e)/;M) 



where there appears the Racah-Wigner 6j symbol of SU (2) and / plays the 
role of the total angular momentum quantum number. Recall that the square 
of the 6j symbol in (|7j) represents the probability that a system prepared in 
the state | ( (a b)^ c) / ; M) will be measured in the state | (a (6 c)e )/ ; M). More- 
over, Racah transforms are the key ingredients to complete the construction 
of the Rotation graph introduced above (note that in this case 

"Rotation" refers to a topological operation on binary tree structures and 
not to rotation operators referred to in ii), see Ref. 2, Appendix A). The 
edge set E = {e} of 0„(V,i?) is a subset of the Cartesian product {V x V) 
selected by the action of elementary j-gates. More precisely, an (undirected) 
arc between two vertices v{b) and v{b') 

e(b,b') = {v{b),v{b')) e (VxV) (8) 

exists if, and only if, the underlying Hilbert spaces are related to each other 
by an elementary unitary operation of the type (jUl). Note also that elements 
in E can be considered as mappings 

{v{b),ni{b)) ^ {v{b'),ni{b')) (9) 

connecting each given decorated vertex to one of its nearest 2(n — 1) vertices 
and thus define a "transport prescription in the horizontal sections" belonging 
to the total space {V x C^"^+^)n of the fiber bundle ((T)). The crucial feature 
that characterises the graph (£>n{V,E) arises from compatibility conditions 
satisfied by 6j symbols appearing in ((7j) . The Racah (triangular) identity, the 
Biedenharn-Elliott (pentagon) identity and the orthogonality conditions for 
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6j symbols |1] ensure that any simple path in 0„(V^, E) with fixed endpoints 
can be freely deformed into any other, providing identical quantum transition 
amplitudes (and of course probabilities) at the kinematical level {cfr. Ref. 2, 
Section 3.1 for more details). 

To complete the description of the structure (V, C^'^"'"^, SU{2y)n we call 
into play M-gates which act on the angular dependence of vectors in 7i;[(b). 
By expliciting such dependence according to 

ni{h) = span{ I b;^,0;JM)„}, (10) 

we write the action of a rotation on a basis vector as 

\b-e\<P'-M'J)^= J2 Di^,{a(3^)\b;e,(j); JM).,, (11) 

where [6, (p) and [6', cj)') are polar angles in the original and rotated coor- 
dinate systems, respectively. D^^j^p (0/^7) are Wigner rotation matrices in 
the JM representation expressed in terms of Euler angles (a/37) which form 
a group under composition |5 (see also Ref. 2, Section 3.2 and Appendix 
Bl, for a discussion on the possibility of selecting "elementary" M-gates in 
some specific cases). In what follows we shall assume for simplicity that, for 
given n and J, it takes just one unit of time to perform one rotation. The 
shorthand notation SU{2y employed in (jT)) actually refers to the group of 
W-rotations, which in turn can be interpreted as the automorphism group 
of the "fiber" C^*^^^. Since rotations in the JM representation do not alter 
the binary bracketing structure of vectors in computational Hilbert spaces 
we can interpret actions of W-matrices as "transport prescriptions along the 
fiber ". 

The framework outlined above makes it manifest the fact that we can 
switch independently j and M-gates without mixing spin and magnetic quan- 
tum numbers. This feature, which relies on the discreteness of the base space 
V and on the "triviality " of the total space {y X C2-^+^)„ , makes it easy with 
kinematics, although the simulator will exhibit complex behaviors at the dy- 
namical level. 

A pictorial representation of the computational fiber space {V, C^"^+^, 
SU{2Y)n for (n + 1) = 4 incoming spins is sketched in Fig. 1. Each vertex 
is associated, through the mapping (j2)), with a binary coupled Hilbert space 
'Hnih) depicted as a plane. Edges of the graph represent Racah transforms 
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(jni) which actually move each space viewed as a whole into one of its nearest 
2{n — 1) = 4 ones. Inside each Ti.^{b) = C^'^^^ (see (jSI)) we can pick up a 
particular vector (as shown in the drawing) and rotate its components by 
means of the Hermitean conjugate of the matrix D'ljj^, (a/57) introduced in 
f|TT|) for some choice of (a/37). 



6 



Figure 1: The fiber space structure of tlie spin network simulator for (n+1) = 
4 incoming spins. Vertices and edges on the perimeter of the Rotation graph 
fSslVjE) have to be identified through the antipodal map. A couple of 
vertices are "blown up" into planes representing the associated computational 
Hilbert spaces. 

The dynamical behavior of the simulator in working as a quantum cir- 
cuit can be discussed in terms of directed paths in the fiber space structure 
{V, C^"^+^, SU{2y)n- By a directed path V we mean a (time) ordered se- 
quence 

l^^in = I'^O )n ^ |t>l )n ■ ■ ■ ^ Ifs )n ^ ■ • ■ ^ \^L )n = It'out ' (^2) 

where we use the shorthand notation \^s)n for computational states and 
s = 0, 1, 2, . . . , L is the lexicographical labelling of the states along the given 
path with fixed endpoints. L is the length of the path P, which turns out to 
be proportional to the time duration of the computation process L ■ t = T 
in terms of the discrete time unit r. The integer L characterising the par- 
ticular directed path in ()12p equals the number of time-ordered elementary 
operations (computational steps) needed to get |Dout)" froin l^iin)" follow- 
ing the path V. An elementary computational step, represented by an arrow 
in (fT^ . is either a Racah transform or a W-rotation. A circuit-type com- 
putation consists in evaluating the expectation value of the unitary operator 
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ilp associated with the path V, namely 

(Oout I I Oin )n- (13) 

By taking advantage of the possibihty of decomposing ilp uniquely into an 
ordered sequence of elementary gates, (fT3|) becomes 

L-l 

(14) 

s=0 

with L = L(V) for short. The symbol [ J-p denotes the ordered product 
along the path V and each elementary operation is now denoted by Us^s+i 
(s = 0, 1, 2, . . . L(V)) to stress its "one-step " character with respect to com- 
putation. Consequently, each elementary transfer matrix in (jl4|) turns out 
to be associated with a local Hamiltonian operator arising from 

(o,+i I D^)^ = exp{iH„(s,s + l) -r} (15) 

and representing the unitary evolution of the simulator in one unit of its 
intrinsic time variable r. We indicate with the shorthand notation (s, s + 1) 
the dependence of on its variables to make clear the local nature of 
this operator with respect to the computational space (V, C^'^^^, SU{2Y)n- 
When (jl5|) is inserted in (jl4j) . such virtual Hamiltonians generally do not 
commute with each other but nonetheless the whole computational process 
may be identified with a well defined unitary evolution of the simulator, 
driven by a poly-local Hamiltonian, in the internal time interval T = L{V)-t. 

Different types of evolutions can be grouped into "computing classes" 
based on the choice of gates that a particular program [i.e. a collection 
of directed paths) has to employ. A computing class that alternates j and 
M-gates would be the most general one (as explained in Ref. 2, Sect. 4.2). 
However, two particular classes to be illustrated below come out to be quite 
interesting in view of the possibility of relating them to other models for 
quantum computation. 

An M-computing class contains programs which employ only (finite se- 
quences of) M-gates in their associated directed paths and it is not difficult 
to realize that such kind of computation, when apphed to N i-spins, repro- 
duces the usual Boolean quantum circuit. 

A j-computing class includes programs which employ only j-gates at each 
computational step. This class is particularly interesting since it shares many 
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features with suitable types of discretized TQFTs, the so-called state sum 
modelsf^jn][H], as we already noticed in Ref. 1 and developed in details in 
Ref. 2, Sect. 5. Now the combinatorial structure of Rotation graphs becomes 
prominent owing to the existence of a 1:1 correspondence between allowed 
elementary operations and the edge set E of the graph 0„(V, E), for any n. 
We denote the unitary operator associated with a program V in this class by 

l^snj ■ |Oin)n — ^ l^^out)" 5 (16) 

where, in a circuit-type setting, |tij^)n is fixed and |0q^^)„ is an accepted 
state. However, in the j-computing class one may address other types of 
problems. For instance, once selected two states, say |0ij2)n and |tiout)".' ^® 
may consider all possible paths V that compute l^out)" result of the 

application of some U^nj to The functional which takes care of such 

multiple choices is a "path sum" (a discretized Feynman's path integral) 
which may be written formally as 

Z[Om> Oout] = Y.^'P ^^OUt I ^3nj I Ojn )n , (17) 

p 

where the summation is over all paths with fixed endpoints and W-p is a 
weight to be assigned to each particular path. The explicit form of such 
functional would contain sums over intermediate spin variables 's of prod- 
ucts of 6j symbols with suitable weights (dimensions of SU (2) irreps labeled 
by kh 's) and phase factors. Note that if we give the same weight, say W-p = 1, 
to each path, then the results on equi-probable amphtudes derived from the 
algebraic identities for 6j symbols ensure that the functional ()17p is a com- 
binatorial invariant, namely it is actually independent of the particular path 
connecting Dj^^ and font- "^^^ situation in this case resembles what happens 
in studying classes of combinatorial invariants (S'?7(2)-state sum function- 
als) for triangulated manifold originally defined in dimension 3 (see Refs. 5,6 
for reviews) and extended (in all dimensions) to manifolds with boundary 
jHJ. Functional like (fT7|) can be actually associated with particular classes 
of extended geometrical objects, namely triangulated surfaces and possibly 
triangulated 3-manifolds. On the other hand, the number of incoming spins 
must be kept fixed in a quantum circuit context, and thus we cannot con- 
clude that such a computing machine is able to simulate TQFTs or gravity 
where the path integrals contain sums over all configurations of fields (up 
to regularization). Other intriguing possibilities open when we insert non 
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trivial weights W-p in (fT7|) . since we may naturally address questions about 
most efficient algorithms and time complexity. For instance we could weigh 
paths with the inverse of their lengths L{V); then the minimum-length path 
(the optimal algorithm) will be dynamically singled out in the path sum. As 
a matter of fact, even such simple example turns out to be highly non trivial 
owing to the combinatorial complexity of (3n{V,E) P, as pointed out in Ref. 
2 (Sect. 4.3 and Appendix A). 

As a final remark we mention the possibility of exploiting the fiber space 
structure (V, C^'^"'"^, SU{2y)n to provide a natural discrete-time implemen- 
tation of holonomic quantum computation ^U]; at least when the Hamilto- 
nians defined in (|15|) exhibit the required degeneracy. 
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